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The ESA's Rosetta spacecraft will arrive at comet 67P/Churyumov-Gerasimenko in 2014. The study of gas and 
dust emission is primary objective of several instruments on the Rosetta spacecraft, including the Microwave 
Instrument for the Rosetta Orbiter (MIRO). We developed a model of dust thermal emission to estimate the 
detectability of dust in the vicinity of the nucleus with MIRO. Our model computes the power received by the 
MIRO antenna in limb viewing as a function of the geometry of the observations and the physical properties of 
the grains. We show that detection in the millimeter and submillimeter channels can be achieved near 
perihelion. 

© 2013 Published by Elsevier Ltd. 


1. Introduction 

Comet 67P/Churyumov-Gerasimenko (67P) was discovered in 
May 1969 after its perihelion at 1.28 AU. Because of repeated close 
encounters with Jupiter in 1840 and 1959, the perihelion distance 
varied from 2.7 AU prior to 1959 to the current value (1.28 AU) and 
the activity of the comet increased. This increase of the activity 
likely facilitated its discovery. Comet 67P was observed during the 
seven perihelion passages since its discovery and an initial 
characterization of the activity (Osip et al„ 1992; Weiler et al„ 
2004; Lara et al., 2011), of the nucleus properties (Lamy et al., 
2007; Tubiana et al., 2008; Kelley et al., 2009) and of the coma 
environment (Hansen et al., 2007; Agarwal et al., 2007, 2010), has 
been drawn. 

The ESA (European Space Agency) Rosetta spacecraft was 
launched on 2 March 2004, to reach comet 67P in May 2014. 
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The probe will be inserted into an orbit around the nucleus at a 
heliocentric distance r h slightly in excess of 3 AU. The lander Philae 
will be deposited on the nucleus in November 2014 (r/,«3 AU). The 
orbiter will be maintained in the vicinity of the comet until 
perihelion (r h =1.3AU) or even until r/,«1.8AU post-perihelion 
(December 2015). Thus, nineteen months of uninterrupted, close- 
up observations of the gas and dust coma will be obtained. The 
Rosetta mission is exceptional because for the first time a space- 
craft will place a lander on a comet nucleus and follow the comet 
during its journey towards the Sun. Scientific instruments on the 
lander and orbiter will characterize the evolution of comet gas and 
dust activity during its approach to the Sun. 

The goals of the Rosetta mission are the global characterization 
of the nucleus (morphology, surface composition, internal 
structure) and of the coma (development of the activity, composi- 
tion and physical properties, dust-gas interaction, interaction 
with the solar wind) with both in situ and remote-sensing 
instruments. The increase of the nucleus temperature leads to an 
increase of the production rate accompanied by the development 
of the coma. 

Rosetta carries a wide range of scientific payload (Glassmeier 
et al., 2007), including the Microwave Instrument for the Rosetta 
Orbiter (MIRO) (Guilds et al., 2007) to characterize the nucleus, gas 
and dust properties of comet 67P. 
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In this paper we present a model of dust thermal emission to 
investigate the detectability of dust with the MIRO instrument as a 
function of the geometry of the observations and of the physical 
properties of the grains. 

2. MIRO instrument 

The MIRO instrument (Guilds et al., 2007) is a consortium 
instrument, built and operated at the Jet Propulsion Laboratory, 
with hardware contributions from the Max-Planck-Institute for Solar 
System Research, and the Paris Observatory. Basic elements of the 
instrument are a 30-cm diameter, offset parabolic reflector telescope, 
and two heterodyne receivers operating at millimeter (190 GHz- 
1 .6 mm) and submillimeter (562 GHzssO.5 mm) wavelengths. The half 
power beam widths (HPBW) of the MIRO beams are 23.8 + 1.5 arcmin 
and 7.5 + 0.25 arcmin at millimeter and submillimeter wavelengths, 
respectively. For a typical comet-Rosetta spacecraft (S/C) distance 
A = 20 km, the spatial resolutions at the nucleus are 
0.138 and 0.044 km at millimeter and submillimeter wavelengths, 
respectively. 

The primary goals of the MIRO instrument are to measure the 
sub-surface temperatures of the nucleus, the gas production rate, 
the relative abundances, the velocity and the excitation tempera- 
ture of gas species, along with their spatial and temporal varia- 
bility. In particular, MIRO has a very high sensitivity to measure 
direct outgassing of three volatile species from the nucleus: CO, 
CH 3 OH, NH 3 , and three oxygen-related isotopologues of water: 
H 2 16 0, H 2 17 0 and H 2 18 0 to obtain the fundamental isotope ratios 
17 0/ 16 0 and 18 0/ 16 0. We refer the reader to Guilds et al. (2007) for 
a detailed presentation of the MIRO instrument. 

The millimeter receiver is configured with a single, broad-band 
continuum detector. The submillimeter radiometer is configured 
with both a broadband continuum detector and a very-high 
spectral resolution (44 kHz) spectrometer. Both continuum chan- 
nels operate in total power mode. The 1 to 5-s Allan deviation 
noise is estimated to be, respectively, 0.11 I< and 0.31 K at milli- 
meter and submillimeter wavelengths (Guilds et al., 2007). These 
1 -a receiver noises will allow us to estimate the detectability of 
dust thermal emission with MIRO. The signal-to-noise ratio of the 
measurement must be > 3c for a statistically reliable detection. 

3. Model of dust thermal emission 

Observations at submillimeter and millimeter wavelengths are 
most sensitive to large particles in the cometary grain size distribution 
and complement optical and infrared observations, which probe 
micrometer-sized grains. By measuring the radiation from large 
particles, we can obtain information on the dust production rate. 
Jewitt and Luu (1990) showed that a substantial fraction of the total 
dust mass in comets is contained within the largest grain sizes. 

Our model computes the thermal emission of the dust coma 
measured in the MIRO beam, expressed in antenna temperature 
scale, which can be compared to the sensitivity of the continuum 
receivers. We assume that the dust coma contains grains with radii 
in the range from a min = 0.1 pm to the maximum liftable size a max . 
The antenna temperature T A is obtained from numerical integra- 
tion along lines of sight in the field of view 

A 2 1 

Ta = rrr-— §B d (d,cp,p)P N (0,<p) dn, (1) 

2K(S £2 a Qm 

where k B is the Boltzmann constant, A is the wavelength, and v is the 
frequency. Pn( 8, cp) is the normalized MIRO beam pattern (i.e., equal to 
1 at maximum) approximated by a Gaussian. Q is the main beam solid 
angle and Q A is the beam solid angle equal to 5.4 x 10 -6 and 
5.4 x IQ 5 sr at 4 = 0.5 and 1.6 mm, respectively. B d (8,cp,p) is the 


spectral brightness (in units of [W m -2 sr -1 Hz -1 ]) of the dust grains. 
We define the spherical coordinate system (z, 6, tp) centered at the 
Rosetta S/C, where z is the radial distance from the S/C along the ray 
defined by {6, cp), and 8 = 0 corresponds to the center of the beam 
(Pn= 1 )■ 

The spectral brightness B d (8,cp,v) is obtained by numerical 
integration along z and grain size a: 

pOO p Qmax 

B d (8,cp,v) = / / F d (a,z,8,cp,p) da dz, (2) 

Jo Ja min 

where F d (a,z,8,cp,v) describes the contribution of particles of size 
a to the emission. In the case of an isotropic coma, F d (a,z,8,cp,v) is 
only dependent on the distance r of the point (z, 4, cp) to the center 
of the nucleus. For thermal emission, the term F d (a, z, 6, cp, p) can be 
written as 

F d (a,z, 8, cp,p) = n d (a,r)xa 2 Q abs (a,p)B,(T d ), (3) 

where n d (a , r) is the density distribution of grains of size a, T d is the 
equilibrium temperature of the grains, and B„(T d ) is the Planck 
function. T d can be assumed in first approximation to be indepen- 
dent of the grain size and taken equal to the equilibrium 
temperature (Jewitt and Luu, 1990). For example we obtained 
T d =244 K and 149 I< at the heliocentric distances r h = 1.3 AU and 
3.5 AU, respectively. The absorption efficiency Qa bs is a function of 
grain size, shape, composition and porosity of the grain. 

We used the Mie theory (Bohren and Huffman, 1983; Van de 
Hulst, 1957), which is the most common technique used through- 
out the literature (Lisse et al., 1998; Hanner et al., 1996, 1994; 
Harker et al., 2002), to compute Qa bs for spherical grains of given 
radius and optical constants m i = n-ik, where n and k are the 
wavelength-dependent indices of refraction. Three grain composi- 
tions are considered: organics, olivine, and two-layer grains. The 
two-layer composition is based on the structure of cometary 
grains outlined by Greenberg and Hage (1990). The basic building 
block consists of a silicate core surrounded by an organic compo- 
nent. The indices of refraction that we used to calculate Q^s for the 
amorphous olivine and organics are from Pollack et al. (1994). 

The Maxwell Garnett effective medium theory was used to 
calculate the effective refractive index for the two-layer grains 
(Greenberg and Hage, 1990). We considered two-layer grains with 
a fractional mass of the mantle component equal to 50%. We adopt 
a density of p 0rg = 1.5 g cm -3 for organics grains and of p ol = 3.5 
g cm -3 for olivine grains (Pollack et al., 1994). 

In order to simplify the analysis we assume that grains of the 
same size have the same mass m d (a) irrespective of their actual 
composition. We assume that the density of grains is p = 1 g cm -3 
which corresponds to a porosity of rs 33%, 71% and 60% for 
organics, olivine and two-layer grains, respectively. We used a 
form of the Maxwell Garnett formula, which is described by Hage 
and Greenberg (1990), to calculate the refractive indices for porous 
grains. The indices of refraction at the wavelengths A = 0.5 mm 
and 4= 1.6 mm used in this study are given in Table 1. 

The density distribution of dust grains in the size range 
(a, a + da), which is a term of Eq. (3), is derived from the continuity 


Table 1 

Optical constants. 


Composition 

>1 = 0.5 mm 

A = 1.6 mm 

Reference 


n 

k 

n 

k 


Olivine 

2.08 

4.00 x 10- 2 

2.08 

3.00 xlO- 2 

Pollack et al. (1994) 

Organics 

2.28 

1.90 xlO- 2 

2.28 

2.80 xlO- 3 

Pollack et al. (1994) 

Two-layer 

2.19 

2.82 x 10~ 2 

2.19 

1.46 xl0~ 2 

This work 
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Table 2 

The nucleus temperature at the sub-solar point and the maximum liftable size for different heliocentric distances and gas production rates. 


r„ (AU) 

T„uc (K) 

(1™°* (kgs- 1 ) 

Q™* (mols- 1 ) 

amax (mm) 

Qf (kgs- 1 ) 

Q™ (mols- 1 ) 

amax (mm) 

1.3 

353 

3.81 x 10 2 

1.13 xlO 28 

16.5 

1.29 xlO 2 

4.14 xlO 27 

5.8 

2.0 

285 

9.10 xlO 1 

2.50 xlO 27 

3.4 

1.58 x 10 1 

o 

X 

o 

LD 

0.6 

3.0 

232 

2.41 xlO 1 

5.40 xlO 26 

0.7 

2.09 

3.70 xlO 25 

5.7 xlO- 2 

3.5 

215 

1.73 xlO 1 

3.80 xlO 26 

0.5 

1.33 

2.10 x 10 25 

3.3 xlO- 2 


equation of the dust flow: 


nd(a, r) da = 


t n(a) da 
4 nr 2 v d {a, r) 


(4) 


where n(a)o<:a~ l, is the size distribution, v d is the grain velocity 
along the radial distance from the center of the nucleus r, Q/ is the 
total dust production rate (in units of [s 1 ]). From the analysis of 
the dust trail of comet 67P the power index q was estimated to be 
3.5 up to 4.1 (Agarwal et al., 2010; Fulle et al., 2010; Ishiguro, 2008; 
Kelley et al., 2008, 2009). For comet 9P/Tempel 1, Gicquel et al. 
(2012) inferred 2.8<q<3.1 from observations with Spitzer under- 
taken before the Deep Impact experiment. We will investigate size 
indexes from 2.5 to 4.5. 

In the following paragraphs, we will describe the parameters to 
compute the dust grain velocity. In our model we assume that the 
distributions of gas density p g and gas velocity v g along the radial 
distance from the center of the nucleus r are described by 
adiabatic spherical expansion of gas from a source of gas where 
the gas velocity on its surface is equal to the local sound velocity. 
This expansion is specified by the radius of the source r n , the gas 
production rate Qg, the surface temperature of the nucleus T n and 
the gas properties: the molecular mass m g and specific heat ratio y 
(Lamb, 1932; Gombosi, 1994). Adopted values for m g and y 
correspond to mixtures of H 2 0, C0 2 and CO gases with relative 
proportions varying with heliocentric distance (Biele et al., 2012). 

We assume that the presence of dust has no impact on the gas 
and we also neglect collisions between dust grains. The motion of 
the grain is governed only by the gas drag force F A and nucleus 
gravity force F c . All other forces (i.e., the solar radiation pressure, 
the solar gravity, etc.) are assumed negligibly small in comparison 
with F a and F c within the domain of simulation. In such condi- 
tions, dust grains with mass m d move radially and variation of the 
dust grain velocity v d along the radius is computed by numerical 
integration of the equations of motion: 



Fig. 1. Grain velocity as a function of distance to nucleus center. The different 
curves correspond to sizes a = 0.1 pm (solid line), a = 10 pm (dashed line) and a max 
(dashed-dotted line). Red and black curves correspond to r ( = 1.3AU and 
0 , = 3.5 AU, respectively. We used the maximum gas production rates (with the 
corresponding maximum liftable sizes) given in Table 2. (For interpretation of the 
references to color in this figure legend, the reader is referred to the web version of 
this article.) 

The maximum ejectable grain size a max is obtained from 
equality of the drag force and the nucleus attraction at the surface 
(Crifo and Rodionov, 1997). Therefore, for sonic expansion of the 
gas and zero initial velocity of the dust grains on the surface and 
taking into account that C D «4 on the surface 


_ 3Qg \J ym g k B T„ 
RnGpM 


dr 

dt 


= v d . 



= F A + F C 


(5) 


with initial condition v d =0 at the nucleus surface. 
The nucleus gravitational force F c is defined as 


F c = -G 


M n m d 

r 2 


( 6 ) 


where G=6.67 x 10 -11 m 3 kg -1 s~ 2 is the gravitational constant 
and M n = 1.3 x 10 13 kg is the mass of the nucleus computed using 
the estimated 67P nucleus radius of r„=2.04 + 0.11 km (Kelley 
et al., 2009) and a bulk density of 370 kg m~ 3 (Lamy et al., 2007). 

The drag force F A is defined as 


F A =yv g -v d ) 2 p g 7ia 2 C D 


(7) 


where C D is the drag coefficient of grains. As the gas production 
rate of comet 67P does not exceed 10 28 s _1 (Schleicher, 2006), the 
minimal mean free path of the molecules in the coma is of the 
order of meters. Since the size of grains is <c 1 m, the gas flow over 
the grains may be considered as free molecular. Therefore, for the 
determination of the drag coefficient C D we use its free molecular 
expression (Crifo et al., 2005). 


To compute the surface temperature of the nucleus T„, we 
supposed that the surface of the nucleus is mainly heated by the 
Sun. The corresponding surface and sub-surface temperatures will 
imply sublimation effects and gas production that will drive the 
velocity of dust. To model the surface temperature of comet 67P, 
we assume here that heat propagates through the regolith by 
thermal conduction only. We neglect the sublimation processes for 
simplicity as we mainly need an estimation of the surface 
temperature at a given heliocentric distance. The diurnal insola- 
tion cycle is computed according to the spin rate derived by Lamy 
et al. (2007), and heat conduction is computed using a 1-D Crank- 
Nicolson scheme over 15 thermal skin depths S th . We refer the 
reader to Leyrat et al. (2011) for a full description of the thermal 
model. The nucleus is modeled as a sphere of radius r„ and the 
temperature is estimated at the equator of the body. We assume a 
bolometric Bond albedo equal to 0.014 (Groussin et al., 2007) and a 
thermal inertia r= 10 J m -2 K 1 s 1/2 similar to the value deter- 
mined for comet 9P/Tempel 1 (Groussin et al., 2007). 

We summarize the nucleus temperature at the sub-solar point 
and the maximum liftable size for different heliocentric distances 
and gas production rates in Table 2. The maximum (Q™ x ) and 
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Table 3 

Dust opacities at 2 = 0.5 mm and 2 = 1.6 mm. 


Compound 

Size index 
(9') 

Size range 
(mm) 

Porosity 

(%) 

*0.5 mm 

(m 2 kg~’) 

* 1.6 mm 

(m 2 kg-') 

Olivine 

2.0 

10- 4 -16.5 

71 

8.0 

X 10- 2 

3.8 xlO- 2 

Olivine 

2.5 

10~ 4 -16.5 

71 

9.6 

x itr 2 

4.0 x 10~ 2 

Olivine 

3.0 

10~ 4 -16.5 

71 

1.2 

x to - 1 

4.1 x 10~ 2 

Olivine 

3.5 

O 

1 

* 

o 

in 

71 

1.5 

X 10-’ 

3.7 xlO- 2 

Organics 

2.0 

1CP 4 -16.5 

33 

8.6 

x nr 2 

1.5 xl0~ 2 

Organics 

2.5 

10~ 4 -16.5 

33 

1.1 

x io-' 

1.6 xl0~ 2 

Organics 

3.0 

10~ 4 -16.5 

33 

1.7 

x io-' 

1.6 xl0~ 2 

Organics 

3.5 

O 

1 

* 

O) 

In 

33 

1.9 

xio-' 

1.1 xl0~ 2 

Two-layer 

2.0 

10~ 4 -16.5 

60 

7.6 

x 10~ 2 

2.3 x 10~ 2 

Two-layer 

2.5 

O 

1 

* 

in 

60 

9.2 

xio - 2 

2.4 x 10- 2 

Two-layer 

3.0 

o 

I 

O) 

in 

60 

1.2 

xio-’ 

2.4 x 10~ 2 

Two-layer 

3.5 

10~ 4 -16.5 

60 

1.3 

x IQ -1 

1.9 xl0~ 2 


minimum (Q™ n ) gas production rates correspond to the expected 
domain of the gaseous activity of comet 67P (Biele et al., 2012). 

We show in Fig. 1 the velocity of the grain as a function of the 
radial distance from the nucleus for various grain sizes and 
heliocentric distances, using the maximum gas production rates 
given in Table 2. At a distance of about 10 nuclear radii («20 km), 
the dust reaches its terminal velocity. 

For comparison with previous studies of dust thermal emission 
at submillimeter wavelengths from ground based observations 
(Jewitt and Luu, 1992; Boissier et al., 2012), we have computed the 
dust opacities (in units of [m 2 kg -1 ]). In first approximation, 
assuming that the dust velocities scale according to v d (a)<xa -05 , 
the dust opacity can be derived as 

_ na 2 Q abs (a, v)n'(a) da 

KX /“"“(4;r/3)^a 3 n'(a) da * 

with n'(a)oca- q and q' = q-0.5. 

Dust opacities at X = 0.5 mm and 1.6 mm are given in Table 3 
for the maximum size expected at perihelion (case Q™ ax ) and a 
range of size indexes. Values are comparable to those determined 
by Jewitt and Luu (1992) and Boissier et al. (2012). 

4. Results of model simulations 

To prepare the planning of the observations, we computed the 
dust emission in the M1RO beam (expressed in antenna tempera- 
ture scale T a ) as a function of the heliocentric distance, for 
different assumptions on the physical properties of the grains 
and dust coma (composition, size distribution). Calculations were 
made in limb viewing geometry, since in nadir sounding the 
measured signal will be dominated by the nucleus thermal 
emission. We expect T A to be higher at perihelion because the 
comet activity is more important. We also expect higher signals for 
low size indexes q since large grains are more efficient radiators at 
millimeter and submillimeter wavelengths. 

The dust thermal emission expected to be detected with the 
MIRO instrument was computed assuming a dust production rate 
Qd equal to the gas production rate Qg (in mass) for all heliocentric 
distances. Fulle et al. (2010) derived a dust mass loss rate of 
~400 kg s _1 at perihelion, a value which is close to the expected 
maximum value of the gas production rate (Table 2). Their dust 
production rate at 3 AU pre-perihelion is also close to Q™ 0 * at 3 AU. 
On the other hand, Fink and Rubin (2012) estimated the dust-to- 
gas ratio for 67P to 0.03-0.5 at perihelion, with a reasonably likely 
value of around 0.1. A dust-to-gas ratio of the order of 1 is 
more consistent with values measured in short-period comets 


(e.g, Colangeli et al., 1998). Calculations were performed for 
heliocentric distances from 1.3 to 3.5 AU. For each distance, we 
considered the two production rates Q™ n and Q™ 1 “ corresponding 
to the expected range of activity of 67P (Table 2). 

4.3. Variation with the geometry of observations 

In this section, we present the evolution of the signal as a 
function of heliocentric distance and limb distance. We consider 
two-layer grains and a size distribution with size index q= 3.5. The 
comet-S/C distance is 4 = 20 km. The signal is computed for a limb 
viewing geometry, and limb distances Su mb between 0.5 km and 
5 km, where 5 Umb is the distance above the surface. Although, the 
thermal emission from the dust increases as the boresight distance 
to the limb decreases, the contribution of the nucleus in the 
sidelobes also increases. The far side lobes of the antenna are not 
well characterized at this time, however it is known that a 
Gaussian beam shape falls off too sharply at large angles from 
the boresight direction. We estimate that the contribution of the 
nucleus in the side lobes ranges from 1 I< to 0.1 I< assuming that 
the MIRO beam can be represented by a Gaussian sitting on top of 
a pedestal whose amplitude is -60 dB to -70 dB below the peak 
gain of the Gaussian. We believe that a conservative estimate is 
that the power pattern has fallen by - 70 dB or more at 0.5 km 
from the nucleus. At the - 70 dB level, the nucleus adds 0.1 K to 
the antenna temperature which is approximately 10% of the 
thermal emission from the dust (assumed to be 1 K). 

The variation of the signal as a function of heliocentric distance 
is shown in Fig. 2 for the two MIRO continuum channels and a 
limb distance Su mb = 1 km. The antenna temperature decreases 
with increasing heliocentric distance, which can be explained by 
the decrease of the dust production rate and the dust temperature 
with increasing heliocentric distance. The antenna temperature at 
X = 0.5 mm is ~3 times higher than at 1.6 mm. 

The variation of the antenna temperature as a function of limb 
distance is shown in Fig. 3 for The antenna temperature 

decreases with increasing limb distance. 



Fig. 2. Antenna temperature for Q™* (black) and Q™" (red) at 2 = 0.5 mm (solid 
line) and 2 = 1.6 mm (dashed-line) as a function of the heliocentric distance. We 
used two-layer grains with a size index q= 3.5. The comet-S/C distance is 
3 = 20 km and the limb distance is Sn mb = 1 km. The horizontal lines show the 
3 -a detection limits obtained for 1 min on-off cycle (blue) and 1 h on-off cycle 
(green), at 2 = 0.5 mm (solid line) and 2 = 1.6 mm (dashed-line). The integration 
time is 5 s on each position. (For interpretation of the references to color in this 
figure legend, the reader is referred to the web version of this article.) 
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Fig. 3. Antenna temperature for Q™" at 2 = 0.5 mm (filled) and X = 1.6 mm 
(empty) as a function of the limb distance. We used two-layer grains with a size 
index q= 3.5. The comet-S/C distance is A = 20 km. The heliocentric distance is 
r h = 1.3 AU. The horizontal lines show the 3-a detection limits obtained for 1 min 
on-off cycle (blue) and 1 h on-off cycle (green), at X = 0.5 mm (solid line) and 
X = 1.6 mm (dashed-line). The integration time is 5 s on each position. (For 
interpretation of the references to color in this figure legend, the reader is referred 
to the web version of this article.) 

4.2. Variation with the physical properties of dust 

In this section, we investigate the variation of the dust thermal 
emission with the physical and compositional properties of the dust. 

The variation of the antenna temperature as a function of grain 
composition is shown in Fig. 4. Results for organics (composition 
index 1), olivine (composition index 2) and two-layer (composi- 
tion index 3) grains are displayed for a size index q = 3.5. The 
antenna temperature at 2=1.6 mm is strongly sensitive to the 
grain composition, and the ratio T a (A = 0.5 mm)/T A (X = 1.6 mm) 
varies from 2 to 8, with the higher value obtained for the more 
absorbing grains (organics). 

The variation of the antenna temperature as a function of the size 
index is shown in Fig. 5 for two limb distances (1 and 5 km) and two- 
layer grains. The relative contribution of small grains increases with 
the increase of q. The antenna temperature decreases with increasing 
q because at millimeter and submillimeter wavelengths, the signal is 
dominated by the emission from large grains. 

The ratios T a (A = 0.5 mm)/T A (i= 1.6 mm) for two-layer grains 
are, respectively, 2.7, 2.8, 3.1, 3.8 and 5.9 for q = 2.5, 3.0, 3.5, 4.0 and 
4.5. The signal at 1.6 mm decreases more rapidly with increasing q 
than the emission at 0.5 mm, as it probes larger sizes. Detection of 
dust thermal emission in the two MIRO channels would allow to 
constrain the dust parameters for a more accurate determination 
of the dust production rate. 


5. Detectability with MIRO 

The detection of thermal emission from dust using the MIRO 
instrument is far from an easy task. As indicated in Figs. 2 and 5, the 
antenna temperatures at 0.5 mm at perihelion and close to the limb 
for the models considered are < 1 K. The antenna temperatures at 
1.6 mm are about three times less. The received power (i.e., the 
antenna temperature) from the dust is added directly to the receiver 
noise power, which is estimated to have noise temperatures of rj 800 I< 
and «3600 K at millimeter and submillimeter wavelengths, respec- 
tively. To further complicate matters, the dust emission produces a 
continuum spectrum which will be difficult to separate from the 



Index for the grain composition 

Fig. 4. Antenna temperature for solid organics (index 1), olivine (index 2) and two- 
layer grains (index 3) at A = 0.5 mm (filled) and A = 1 .6 mm (empty). The size index 
is q= 3.5. The comet-S/C distance is 4 = 20 km. Calculations are for Q™ ax at 
r h =1.3AU, and a limb distance <5 /imb = l km. The horizontal lines show the 3-tr 
detection limits obtained for 1 min on-off cycle (blue) and 1 h on-off cycle (green), 
at X = 0.5 mm (solid line) and X = 1.6 mm (dashed-line). The integration time is 5 s 
on each position. (For interpretation of the references to color in this figure legend, 
the reader is referred to the web version of this article.) 



Fig. 5. Antenna temperature as a function of the slope of the distribution at 
A = 0.5 mm (filled) and A = 1.6 mm (empty) for two-layer grains. The comet-S/C 
distance is 4 = 20 km. We used Qj 0 * at r h =1.3 AU. The limb distance is <5 /imf) = l km 
(black) and 8n mb = 5 km (red). The horizontal lines show the 3-cr detection limits 
obtained for 1 min on-off cycle (blue) and 1 h on-off cycle (green), at A = 0.5 mm 
(solid line) and A = 1.6 mm (dashed-line). The integration time is 5 s on each 
position. (For interpretation of the references to color in this figure legend, the 
reader is referred to the web version of this article.) 

continuum spectrum of the nucleus in the side lobes of the antenna 
because of the poor signal-to-noise ratio. 

Separating the signal from the dust from the radiometer noise 
requires a position switching strategy. As can be seen in Figs. 2-5, 
the antenna temperature due to the dust decreases rapidly as the 
limb distance of the boresight increases away from the limb. For 
example, in Fig. 3, the antenna temperature at 0.5 mm decreases 
by «0.5 K between boresight pointing positions of 0.5 to 5 km. The 
l-o- receiver noise of the submillimeter receiver is estimated to be 
«0.31 K (0.11 K for the millimeter channel) in a 5 s integration 
time. By position switching between a boresight on the dust and blank 
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sky, and treating the temperature difference between T^(dust) and 
T/i(sl<y) as an independent measurement, it should be possible to stack 
many measurements and to greatly improve the signal-to-noise ratio. 
In a 10 s on-off cycle (5 s on each position), the noise on a differential 
submillimeter measurement will be 0.4 K (0.02 K in an hour). This 
yields a signal-to-noise ratio of 15 based on a 3 -a noise level and a 1 K 
contribution from the dust. Another alternate approach to position 
switching will be to perform repeated continuous scans away 
(or towards) from the limb and perpendicular to the limb. This 
observing mode should be more efficient than position switching, by 
avoiding start, stop, and settling times associated with position 
switching. This would allow a greater signal-to-noise ratio for the 
same integration time mentioned above. For data acquired close to the 
limb edge, it may be possible to subtract the nucleus thermal emission 
from the individual scans or stacked scans. 

In Figs. 2-5, are plotted the 3 -a detection limits obtained for 1 min 
and 1 h on-off cycle for comparison with the expected antenna 
temperatures at 0.5 and 1.6 mm. Detection of dust above the 3-a 
threshold is expected at perihelion for both continuum channels, 
providing the size index q < 4.0, which is in the range of expected 
values for 67P (Fulle et al., 2010). The signal-to-noise ratio is not high 
enough to provide significant constraints on the size distribution, but 
might allow us to exclude extreme cases, such as a very steep size 
distribution or very absorbing material such as pure organic grains. 

A dust-to-gas ratio higher than 1 would be more favourable for 
M1RO detection of dust. In addition, the local gas and dust production 
rates on certain regions of the nucleus surface (e.g., active areas or 
subsolar regions) will be probably higher than those considered in our 
model which assumes isotropic production of dust and gas. However, 
worth is to note that keeping the dust-to-gas ratio equal to 1, the 
antenna temperature is not proportional to the dust or gas production 
rate. We made computations with the parameters of Fig. 2, by 
assuming at r h = 1.3 AU a gas production rate 10 times higher than 
Qg 0 * in Table 2. With respect to results obtained with Q™ 0 *, the 
antenna temperature increases by a factor of 1.6 and 3.5, at 
X = 0.5 mm and 1.6 mm, respectively. This behaviour can be explained 
by the 10 times increase of the maximum liftable size (now 17 cm). 
Most of the dust mass is then in large particles which have a small 
emitting cross-section. 

6. Conclusion 

We modeled the thermal emission of dust grains in comet 67P in 
order to investigate the detectability of the dust coma by the MIRO 
instrument aboard Rosetta. We show that detection in the millimeter 
and submillimeter channels can be achieved near perihelion. MIRO 
observations will provide constraints on the dust production of 
millimeter-sized particles which might complement measurements 
obtained from the Grain Impact Analyser and Dust Accumulator 
(G1ADA) experiment (Colangeli et al., 2007). MIRO data might also 
be combined to infrared measurements of dust thermal emission 
performed with the Visual IR Thermal Imaging Spectrometer (VIRTIS) 
(Coradini et al., 2007) for further insights into dust properties. 
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